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Abstract 

We analyze the behaviour of moduli fields in string effective models between 
the end of inflation and reheating. The effective moduli potential during this era 
is derived for a class of simple models. We argue that this potential significantly 
stabilizes the modulus at its high energy minimum, if some restrictions on modular 
weights are met. Two mechanisms to further stabilize the moduli to their low 
energy minima are discussed explicitly: coinciding minima at a point of enhanced 
symmetry, and the smooth transition from high to low energy minimum by an 
effective mass term C 2 H 2 . For both cases we present explicit examples, and C 2 is 
found to be O(10) at most. In addition, we show that during a smooth transition 
the reduction of the modulus amplitude strongly depends on the shape of the low 
energy potential. 
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1 Introduction 



In view of the current understanding of superstring theory, moduli fields are likely to have 
masses of the gravitino mass scale, which is assumed to be at, or at least not far from, 
the electroweak scale. The existence of such very weakly interacting scalar fields with an 
electroweak-scale mass would cause a cosmologically disastrous problem, which is referred 
to as the (cosmological) moduli problem in the context of string theory ||T|. This problem 
can be regarded as a string version of the notorious Polonyi problem H associated with 
a hidden sector supersymmetry breaking in supergravity. The most severe consequence 
is that late moduli decay after coherent oscillation would totally change the abundances 
of light nuclei successfully predicted by the standard big-bang nucleosynthesis theory. 

A number of ideas to solve the moduli problem (and the Polonyi problem) have been 
proposed over the last decade, some of which have been reviewed in ref. ||. What makes 
the problem particularly serious is the fact that unlike most other unwanted relics the 
energy stored in moduli cannot be diluted away by ordinary inflation. At the end of 
inflation moduli will generically be displaced by 0(1) in Planck units from their low 
energy minimum and finally (once the Hubble constant approaches their mass) start late 
decaying coherent oscillations with the aforementioned disastrous consequences. Two 
types of (cosmological) solutions have been proposed in view of this particular problem. 
The first assumes a second late period of inflation [|J (e. g. at the weak scale) of a few 
e-foldings sufficiently long to dilute away the energy density stored in the moduli. A very 
interesting realization of this mechanism is provided by the idea of thermal inflation ||. 

In this paper we will be concerned with the second type of solutions which attempts 
to stabilize the moduli to their minimum before the coherent oscillations start. A useful 
observation towards such a solution is that all fields receive an effective mass of the order of 
the Hubble parameter H due to the energy density present in the universe JJ. This mass, 
which will dominate in early stages can stabilize moduli very effectively to a minimum. 
Unfortunately, this minimum - the high energy minimum - will generally not coincide 
with the low energy minimum and oscillations reappear as soon as the moduli start to 
roll towards the second. It has, however, been suggested that symmetries might provide 
a reason for the low and high energy minima to coincide J7J. 

In string theory the role of this symmetry can be played by modular invariance with 
the minimum sitting at a point of enhanced symmetry. In fact, a self-dual point is an 
extremum of the potential and the moduli problem can be solved once this extremum is 
arranged to be a minimum in each stage of the cosmological evolution. Though possible in 
principle, this scenario poses a problem when it is realized during the inflationary epoch: 
In models based on low energy effective string actions inflation is generally difficult to 
obtain || and it appears to be problematic to construct explicit examples as long as 
the precise mechanism of inflation is unknown. Not a single consistent example with an 
inflationary behaviour which stabilizes a modulus to a point of enhanced symmetry is 
known to us. 

Recently Linde || has proposed a dynamical damping mechanism to reduce the am- 
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plitude of the moduli which does not assume coinciding high and low energy minima. 
The crucial observation is that an effective moduli mass term % ie g ~ C 2 H 2 /2 with a 
large number C 2 ^> f causes a significant damping of the oscillation amplitude. This 
happens because such a large value C 2 guarantees a 'smooth' transition of the modulus 
from the high to the low energy minimum at the time when H has decreased down to 
the low energy modulus mass rnx and the low energy potential becomes relevant. In 
ref. it has been shown that this damping mechanism works for a quadratic low energy 
potential Vq ~ m^T 2 /2 and leads to a damping of the oscillation amplitude by a factor 
of ~ cxp(— irpC/2) where p described the expansion rate of the universe via H = p/t. 

In this paper, we would like to study the moduli potential in the post inflationary era 
with particular emphasis on the aforementioned possibilities to solve the moduli problem. 
In particular, we will derive the effective scalar potential for a modulus which describes its 
behaviour during the preheating phase after inflation, namely the era when the inflaton 
obeys a damped oscillation until it decays to reheat the universe. It is this era that 
determines the initial condition of the moduli oscillation. Therefore, for a solution of the 
moduli problem based on a symmetry induced minimum this era has to be considered 
necessarily. Moreover, we do not have to deal with consistency problems which arise 
from the unknown mechanism of inflation in string effective models if the stabilization of 
moduli is implemented in the preheating phase. All we have to assume is the existence 
of such a phase with the energy density being dominated by a coherently oscillating field 
which - depending on the type of inflation - might or might not coincide with the inflaton. 
Indeed, we will construct an explicit example which shows the stabilization of a modulus 
at a self-dual point during this era. 

Also the mechanism of ref. |§ mentioned above can be addressed in our framework. 
A solution of the gravitino problem ]nj implies a low reheating temperature ~ 10 9 GeV, 
corresponding to the Hubble parameter ~ 1 GeV, well below the electroweak scale. Since 
the smooth transition between high and low energy minima occurs once H drops to 
H ~ rnx it will typically take place during the preheating era. In our framework, we are 
also able to explicitly compute the value C 2 which was put in by hand in ref. || and it 
can be checked whether it can be really made as large as required. 

To illustrate the above ideas, we will concentrate ourselves on a simple class of models 
consisting of one modulus T with a low energy potential Vo and an oscillating field (in- 
flaton) They will be introduced in the next section where a general analysis of their 
properties during preheating is presented. In particular, we derive an effective modulus 
potential by performing a short time average over one period of the $ oscillation. In 
section 3 we will focus on string motivated models: The low energy potential is assumed 
to be dual and to originate from gaugino condensation [|TT|] and the (modulus dependent) 
Yukawa coupling of $ is chosen in accordance with modular invariance. For these models 
we will first discuss the mechanism of stabilization at the high energy minimum in some 
detail. It turns out that this first stabilization (which is essential for the other mechanisms 
to work) is possible for a restricted range of modular weights of $ only. Then we construct 
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an explicit example with coinciding high and low energy minima at the self-dual point. 
As for Linde's mechanism, we first perform a 'phenomenological' analysis for T-dual low 
energy potentials. We find that the conclusions concerning a suppression of moduli oscil- 
lations can differ substantially from the case of a purely quadratic potential considered 
in ref. ||. After this, information about the possible values C 2 is extracted. It turns 
out that large values are not associated with large but with steep couplings in the theory. 
In this respect the properties of modular functions naturally present in string effective 
models turn out to be very interesting. The results will be summarized and commented 
on in section 4. 



2 The general framework 

In this section we would like to present a class of simple models suitable to describe the 
behaviour of moduli fields during the preheating period. We will investigate their general 
properties and in particular derive an effective moduli potential which incorporates the 
effect of the energy density stored in the coherent oscillations. Explicit examples will be 
given in the next section. 

Our models consist of two chiral superfields, namely a modulus T and the oscillating 
field We do not restrict ourselves to any specific type of inflation which precedes the 
preheating phase. For chaotic inflation |I2] $ typically will be the inflaton itself whereas 



for hybrid inflation |13[ a false vacuum expansion is followed by oscillations of $. 

The model is specified by a Kahler potential and a superpotential of the following 
formQ: 

K = K (T,f) +Z(T,f)|$| 2 



W = W (T) + ±F(T)$ q 



(1) 



Here the Kahler potential K of the modulus and the metric Z of $ are kept general. 
In the examples we will concentrate on string motivated expressions for these quantities. 
Throughout this paper, we use Planck units M ~ 2.4 x 10 18 GeV = 1. We have neglected 
higher $ terms in eq. (ffl) which may dominate the potential during inflation, because we 
focus on the preheating phase where $ < 1. For simplicity, we assume that the $ field 
oscillates around the origin. The superpotential has been split into two parts: a modulus 
potential Wq originating from nonperturbative effects like gaugino condensation and a 
$-term with a modulus dependent coupling F. We allow the power q to take the values 
q = 2, 3 corresponding to a mass term or to a Yukawa coupling for $. Let us denote the 
intrinsic mass scales in Wq by = 0(Wq) and in F by Mq> = 0(F). To be consistent 
with the constraint on the anisotropy of the photon background the mass (coupling) M$ 
should be well below the Planck scale, typically M$ ~ 10 _5 -lCr 6 for chaotic inflation. 
For hybrid inflation, M$ can be taken somewhat smaller. We will assume throughout this 
discussion that the low energy modulus mass ttlt is much smaller, i. e. itlt <C M$. For 



1 Chiral superfields and their scalar components will be denoted by the same symbol. 
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moduli potentials based on gaugino condensation we have my ~ m 3 / 2 ~ 10 15 so that 
this can be safely assumed. 

Generally, the behaviour of the model (|l|) is expected to be complicated. It has to be 
analyzed in the full field space of four real fields and using the complete scalar potential. 
Having the ordinary scenario of reheating in mind one might think about the following 
situation: A rapidly oscillating $ dominates the energy density of the universe and the 
modulus is moving in this given 'background'. Then an effective potential for the modulus 
can be obtained by a short time average over one period of the $ oscillation. Such an 
average procedure would simplify the situation considerably. Let us now analyze under 
which assumptions this picture is really correct. 

As soon as $ 1 the scalar potential can be expanded as 
V ~ Vo + Vx 

Vi = f{T,T)\<S>\ 2q - 2 (2) 



/ = Z~\ K «\F\ 2 



where Vq is the low energy modulus potential originating from Wq. Higher powers in $ 
have been neglected. First we should require that the energy density is dominated by 
i. e. Vo C Vi. Then H 2 ~ V\ and the condition is consistently fulfilled as long as 

H 2 > V ~ m% . (3) 

Since the low energy mass wit is much smaller than the inflationary scale Mq> this condition 
will be fulfilled for a wide range in H. Even for H ~ tjit it may hold if T is close enough 
to its low energy minimum. We assume that this is the case in view of the fact that we are 
concerned with a solution of the cosmological moduli problem where the amplitude of the 
moduli oscillation is either small or damped quickly. Self-consistency of this assumption 
should be checked in each case. From the expressions (0) we can estimate the effective 
masses for T and $ as follows: 



"4, eff = \ f j?H 2 (4) 
H 2 



with 

A short time average is justified if rriT,cS <S ^$,eff which in turn implies 

(7) 

JT 
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Inspection of eq. (|6]) shows that this condition is usually fulfilled once |$| is sufficiently 
small. However there is one important exception: The quantity faff can become very 
large close to zeros F(T) = of the Yukawa coupling. For the following derivation we 
therefore assume that the system - at least for a certain period after |$| became small 
does not evolve towards such a point or, equivalently, that the coupling strength \F\ is 
bounded from below in the relevant part of field space. Later on, we will comment on 
what happens if this condition is violated. 

Oscillation of $ will dominate the energy density of the universe as long as H = p/t ~ 
T$ ~ Tr H where H = a/a is the Hubble parameter, T$ is the decay constant of $ and 
Trh is the reheating temperature. For a radiation like behaviour of H we have p = 1/2, 
whereas a matter like behaviour is characterized by p = 2/3. In view of the gravitino 
problem Trh ~ 10 9 GeV. This implies that the universe can be dominated by $ oscillations 
down to values of H ~ 1 GeV. In particular the start of the coherent modulus oscillations 
in the low energy potential responsible for the moduli problem at H ~ ttit can fall into 
this era. 

Let us now carry out the averaging procedure to obtain an effective potential for the 
modulus in the background of the oscillating field $. As discussed, we assume that the 
conditions (|3|) and ([?[) hold. The system is described by a Lagrange function 



£ = a 3 
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with the Kahler metric K ik = did k K and the fields denoted by <ft l . For simplicity, we have 
assumed a flat Robertson- Walker universe. The conjugate momenta are given by 

vr, = a 3 K fk $* (9) 

and the equations of motion can be written in the form 

j t (K ih p) + SHK&P - (d- k K m )^$% d k V = 0. (10) 
Sometimes it is useful to normalize the second derivative term in eq. fll0|) to get 

f + 3fff + r; fc #' + w = o, (ii) 

with the Kahler connection T defined by = K kl diKjj. Furthermore, the Friedmann 
equation reads 

H 2 = K ik ^$ k + V . (12) 
A straightforward computation under the above assumptions and using eq. (0) leads to 

l(7r*$)~a 3 (l7-(<7-:L)Vi) , (13) 
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where U = K i ^ i ^ k is the kinetic energy. The short time average (which we denote by 
< .. > in the following) of d(ir $$)/ dt approximately vanishes and we have 

<U >~(q-l) <Vx> . (14) 

With the energy density and the pressure p =< U + V >, P =< U — V > we derive an 
equation of state 

P = vp , v = . (15) 

Q 

and the Hubble parameter behaves like 

H=- 2 - 9 (16) 

t ' P 3(l + i/) 3g-3 ' 1 ' 

As expected, for q = 2 (quadratic $ potential) the expansion corresponds to a matter 
dominated universe and for q = 3 (quartic $ potential) it corresponds to a radiation 
dominated universe. Since < V% >= H 2 /q we can carry out the following replacement in 
the EOM for T 

dVl Iff rr2 / 17 X 

< W > - qj ■ 
This leads to the effective EOM for T during preheating 

f + 3Hf + l^r f2 + K ^ d ~w = > ( 18 ) 

with the potential V e g given by 
Kff = V Q + V l = V + ±gH 2 

(19) 

g = In/ = - \nZ + K + \n\F\ 2 . 



For H 3> my the potential ( |I~9"D is dominated by the second term originating from the 
$ oscillations. The modulus will therefore settle down to one of the minima T\ of g. This 
value can be assumed as an initial value for the further development of T which sets in 
once the low energy potential becomes important, i. e. once H approaches tut- In the 
spirit of ref. we can introduce a parameter C 2 

C 2 (T) 1 9T 



qT-Ti 



9t = -^r + K 0T + ^. (20) 
Then the EOM can be rewritten as 



f + 3HT + E^IZf 2 + K-' f | |p + C 2 (T)H 2 (T - T x ) ] = . (21) 



6 



Notice that, in contrast to ref. 0, C 2 depends on T and can - along with the starting 
value Ti - be computed explicitly once a specific model has been chosen. 



Let us now discuss the general consequences of this result before we turn to concrete 
models. Of course it might happen that the system runs into a region of very low coupling 
|F(T)| even before the above stage is reached, i. e. before |$| is sufficiently small. In this 
case which strongly depends on the initial conditions there is no chance to generate a 
substantial damping of the modulus and we are left with the usual moduli problem. 

If however the system evolves according to the above scenario, a damping of the 
modulus oscillation might result. First we would like to analyze the stabilization of T to 
the high energy minimum Ti. Clearly, gr{Ti) = 0, and for Ti to be a stable minimum the 
requirement 

g TT {T x ) > |0zr(Ti)| (22) 

has to be fulfilled. In particular g T f > 0, which might impose severe restrictions on 
the model since g T f does not depend on the superpotential but on the Kahler potential 
only. The stability condition (p2"|) can be used to derive an upper bound on the value of 
C\ = C 2 (T 1 ) which determines how effective the stabilization to Ti can be. From eq. ( pUD 
we have 



C 2 (T ~ Ti) ~ - 
Q 



rji rp 

9tt(Ti) + T _r r 1 9Tr(Ti) 



(23) 



and with condition (22 ) it follows that 

ICfl < - q g T f{Tx) . (24) 

Moreover, the order of \Cf\ will be given by the RHS of eq. (p4|). Notice that again this 
value depends on the Kahler potential of the theory only. The implications of eqs. ( 22[) 



and ( |2~3| ) for string motivated models and the stabilization process will be studied in more 
detail in the next section. 

For the moment we assume that T has been driven to Ti with a certain precision and 
starts its further evolution from this point. If Ti coincides with the low energy minimum 
To the stabilization is definite and moduli oscillations are suppressed. 

Suppose that Ti and To do not coincide. Then, once H ~ vtlt the modulus will feel 
the low energy potential and will start oscillations around the true minimum T . The 
moduli problem simply reappears. In ref. || it has been shown that assuming a quadratic 
low energy potential these oscillations can be suppressed substantially for large values of 
C 2 by a factor of order exp(— npC/2). This mechanism comes into operation as soon as 
H drops to H ~ my. Therefore, it may play a role for the models considered here if 
Trh ~ vrij. 2 so that domination of the $ oscillations is maintained sufficiently long. Such 
a low value of the reheating temperature is also required in order to solve the gravitino 
problem. 
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Large values of C 2 cannot be obtained by simply introducing a large coupling in the 
theory As eq. (|20D shows a coupling constant in front of F or Z simply drops out. Instead, 
what one needs is a very steep coupling such that Ft/F becomes large. Since generally C 2 
is proportional to V'/V the validity of this statement is not restricted to our particular 
class of models. In the next section we will discuss how the conclusions of ref. M change 
for T-dual low energy potentials and to what extent large values of C 2 can be obtained 
in string motivated models. 

Finally, the modulus might run into a region with vanishing Yukawa coupling |-F(T)| 
even after our effective potential description becomes valid. As can be seen from eq. ( |19|) 
these points show up in the effective potential as singularities V e s — > — oo. If such zeros 
exist the system might therefore easily be attracted to them. To analyze the qualitative 
behaviour close to such a point T 2 let us expand F ~ F X (T — T 2 ). Then, according to the 
condition (0) our averaging procedure breaks down once 1/C 2 ~ \T — T 2 \ ~ |$| and T 
will start to oscillate around T 2 with an initial amplitude of |$|. Since |$| can be quite 
small when this happens a large suppression of the T amplitude might result. Again, this 
guarantees a (at least partial) solution of the moduli problem only if T 2 coincides with 
the low energy minimum T . 

For comparison we would like to briefly mention the requirements on the amplitude 
of the coherent moduli oscillation posed by the nucleosynthesis constraint. The energy 
density pr of these oscillations at the time of the decay of T divided by the photon density 
n 7 is given by 

— - TrhT 2 , (25) 

being the initial oscillation amplitude and assuming a matter dominated preheating 
phase. 

In order not to dissociate light nuclei this quantity should be less than typically 
to 10~ 12 GeV (for ordinary moduli masses). For a reheating temperature of Trh = 10 9 
GeV this is satisfied, if T- m ~ 10~ 9 to 10~ n in Planck units. 

In the case of a quadratic low energy potential as used in |J one therefore needs values 
of C 2 £ 10 2 - 10 3 . 



3 Applications and examples 

In this section we would like to discuss the various possibilities mentioned above more 
explicitly. We will concentrate on string motivated models specified by 

K = -w\n(T + T) , S= * . (26) 

Usually, for a modulus we have w = 3. The modular weight of $ is denoted by n, 
i. e. under T — > 1/T duality $ transforms as $ — > $/T n . To have a duality invariant 



S 



theory the modular weight m of the Yukawa coupling F should be given by m = w 
Application of eq. (|19D leads to an effective potential 



g = {n-w) ln(T + T) + In |F| 2 , 

and the quantity which determines the value of C 2 is given by 

n — w Ft 

Qt = = H • 

y T + T F 



qn. 



(27) 



(2? 



The potential g is modular invariant if and only if $ has a vanishing weight n = 0. Since 
$ has been 'integrated' out to arrive at the above result this property could be expected. 

The explicit form of Yukawa couplings in string theory has been computed in a number 
of publications [14], I5| . Generally, it has been found to be consistent with modular 
invariance. As an illustrative class we will consider couplings of the form 



F(T)$ 9 = M$ 



r] 2m (T) 



(29) 



with the Dedekind function 77 and the totally modular invariant function j. By symmetry 
requirements F could contain an arbitrary function of j(T). Here we study only the 
subclass ( pj| ) with an arbitrary power I of j. Yukawa couplings containing j-functions 
show up for example in Calabi-Yau compactifications |l6j . In orbifold models the modular 
weight n of $ has been found to be in the range flT7 



10 > n > -4 



(30) 



Furthermore, we concentrate on a specific form of the low energy moduli potential, 
which is related to the potential one expects from the process of T-dual gaugino con- 



densation [18[]. Assuming gaugino condensation in a model consisting of dilaton S and 
modulus T in the usual way, i.e. 



K 
W 



-ln(S + S) -31n(T + T) 
n(S)/ V (Tf , 



one arrives at the scalar potential 
1 



V Q (S,T) 



TMT) 



121 



A(S) + B(S) 



2 W(T) +1 
V (T) 



(31) 
(32) 



(33) 



where Tr = T + T. The functions A(S) and B(S) are given in terms of Q(S). Usually 
A(S) = —B(S) {B(S) > always) at the minimum (supersymmetry is not broken by the 
dilaton, but by the modulus) and this will give the famous value T w 1.23. The shape of 
this potential is shown in fig. 1. 



9 



Fig 1: Gaugino condensation potential with minimum at T = 1.23 



If B(S) = 0, which corresponds to unbroken supersymmetry in the T-direction the 
minimum will be at T = p (the other self-dual point with p = exp(z7r/3)). This implies 
that supersymmetry is unbroken by this condensation process since for the factorizing 
Ansatz (^) supersymmetry cannot be broken by the dilaton. The potential looks like in 
fig. 2. 
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Fz^ Gaugino condensation potential with minimum at T = p 



A stabilization of the dilaton is usually difficult to achieve and we will not address 
this question here. Instead, to keep our models as simple as possible we will consider only 
the T-dependence of the potential. Consequently, we assume that the dilaton VEV is 
fixed by some unknown mechanism and A(S), B(S) are thus determined. 
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Of course it is well known that the potentials discussed here generically have a negative 
cosmological constant. Sometimes one can shift the vacuum energy to zero by adjusting 
the value of an additional constant term in the superpotential. We simply assume that 
the potentials we study are shifted in such a way that the cosmological constant vanishes. 



3.1 Stabilization at the high energy minimum 

Now that we have specified our models let us apply the general results of the previous 
section. First we study the stabilization to the high energy minimum Ti, which is essential 
for any further suppression of moduli oscillations. As we have seen T\ is a stable minimum 
of the high energy potential provided eq. (|22|) holds. Since 

w — n 

9tt 



(T + T) 2 ' 

this implies as a necessary (but not sufficient) condition that 

n < w = 3 , (34) 

which constitutes a significant restriction on the modular weight n of <£>. If this bound is 
violated, no stabilization at the high scale and consequently no further stabilization can 
take place. Furthermore, from eq. (|24]) we have 

2 \w-n\ 

Cl -gm + fO 2 ' (35) 

Given the bound (|30|) on the modular weight and the fact that the starting value Ti will 
be typically 0(1) this restricts C\ to a few at most. In view of this constraint it seems 
unlikely that the moduli problem can be solved completely by the 'C 2 -mechanism' of 
ref. || though a larger value of C 2 (T) at T ^ T\ cannot be excluded. We will come back 
to this point later on. 

The stabilization of T to 7\ is governed by the EOM Q (cf. eq. (|21"|)) 

(T + T ) 2 

f + 3HT + w ClH 2 {T - Ti) ~ . (36) 

Note that we have omitted the low energy potential V$ which is unimportant at this early 
stage. The solutions of eq. (^) can be written as (T — Ti) ~ t~ a with 




Damping of the T oscillations 
AT f fU\ a fm T 



AT; \tfj \M, 



^ r 't ri) q ■ (37) 



(38) 



2 The term proportional to T 2 in eq. (|T]) can be consistently neglected. 
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and typical maximal damping ratios are 10 5 and 10 2 5 for p = 2/3 and p — 1/2 respec- 
tively (for M # ~ 10" 5 and rriT ~ 10 15 ). Though this suppression itself is definitely not 
sufficient to solve the moduli problem for AT; ~ 0(1) it significantly improves a preex- 
isting suppression AT; <C 1 originating from the inflationary epoch. Of course one has to 
make sure that no destabilization occurs in the following evolution which is the problem 
we are going to address now. 

3.2 High and low energy minima coincide 

Models in which the minima of the high energy potential and the low energy potential 
coincide, are very attractive from the viewpoint of the cosmological moduli problem. In 
this case the modulus will be fixed at this minimum at early times, when the curvature 
at the minimum is large as discussed in the previous subsection. During the transition to 
the low energy potential it will safely stay at that point. 

In principle, this situation can be achieved by carefully tuning parameters in the 
potential. Such a choice would be neither natural nor stable under radiative corrections 
and is clearly unreasonable. Therefore, coincidence of low and high energy minima should 
be guaranteed by a symmetry. In a string framework duality is a very natural candidate 
for this symmetry which has been proposed in ref. |7J . 

For our models, this can be realized for n = only since otherwise duality is broken 
by the effective potential. This fixes the modular weight of the Yukawa coupling F to 
m = w = 3. Then, by modular invariance low and high energy potential both possess 
extrema at 1 and p = exp(i7r/3). For a successful example, one of these points should be 
arranged to be a minimum for both contributions to the potential. 

For T\ = p we can make a more detailed prediction of the stabilization at the high 
energy minimum. The second term in eq. ([23|) vanishes at 7\ = p and the value of C 2 at 
this high energy minimum can be explicitly computed to be 



corresponding to a typical damping of 10 3,3 and 10 2,5 for p = 2/3 and p = 1/2, respec- 
tively. 




1 



(40) 




(41) 



12 



We start with an example for a minimum at p. A simple Yukawa coupling with weight 
m = 3 is given by 



M, 







rf(T) 



(42) 



The resulting effective potential g has a minimum at T\ = p. As discussed above, the low 
energy minimum of Vq is at p if supersymmetry is unbroken (because then B(S) = 0). 

A model where the minima coincide at the other self dual point T = 1 can be con- 
structed as well. A Yukawa-coupling 



M$j l (T) 
V 6 (T) 



-0.068 > I > -0.016 



(43) 



will give the required minimum of g at T = 1. For a superpotential of the form Wq = 
£l(S)/ \rf{T)j(T)) instead of Wq = fl(S)/r] 6 (T) the low energy minimum in T-direction 
will be also at T = 1. The low energy potential of this model looks like in fig. 3. 




Fig 3: Potential for W = tt(S)/(r) 6 (T)j(T)) 



Although explicit orbifold calculations do not lead to j-functions |L5L this superpo- 
tential is perfectly consistent with duality and might originate from Calabi-Yau compact- 



ifications, where similar moduli dependencies show up ||T6|. 

Clearly, stabilization of a modulus at one of its self-dual points implies unbroken 
supersymmetry in this direction. Thus, if supersymmetry is broken by a modulus at all 
at least this modulus cannot be stabilized in this way. 



3.3 High and low energy minima do not coincide 

The situation described above may not be expected for all the (numerous) moduli of a 
superstring theory. Most probably not for each modulus the minima of the low and high 
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energy potentials will coincide. We should therefore study this case and analyze, whether 
the mechanism described in will sufficiently damp the coherent oscillations during the 
transition of the modulus from T\ to To (its final minimum). First we will generalize the 
approach of |§ to a more realistic low energy potential, but keep the high energy potential 
quadratic. As we will see, this will already change the results significantly. 

3.3.1 High energy minimum coincides with maximum of low energy potential 

A smooth transition does not take place at all if the high energy minimum coincides with 
a low energy maximum. 

At this point we have V^{T X ) = V{(Tx) = and > 0, V" < 0. It is clear that at 
all times V^iTi) = and that at some point of time V^j(Xi) = 0. Therefore the position 
of the minimum of V e s will not move smoothly from the initial to the final low energy 
minimum with increasing time. Instead, the modulus will be released from its minimum at 
Ti at the time when V^(Ti) = 0. Thus no reduction of the coherent oscillation amplitude 
can result, regardless of how strong the high energy potential is (i.e. regardless of the 
value of C 2 ). The mechanism of 11 cannot work in this case. 

One can show that this situation occurs not only if the early minimum and the late 
maximum coincide, but also when they are close to each other. The maximal separation 
of the two extrema depends of course on the shape of the potentials, but is usually rather 
small. 

Nevertheless this behaviour gives serious restrictions for our gaugino condensation 
inspired model, since the low energy potential has a maximum (in the real direction) at 
the self-dual point T — 1, if supersymmetry is broken in the T-direction (as is usually 
the case). Therefore high energy potentials whose minimum is at T = 1 will not succeed 
in damping the coherent oscillations of the modulus. 

Fig. 4 shows the trajectory of the modulus, where V\ = C 2 H 2 (T — 1.001) 2 /2 for C — 10 
(a rather large value). It is clear that there is no damping at all, when compared to the 
trajectory with 7\ = 1.7 and C — 10, which is also given in fig. 4. 

3.3.2 Smooth transition between high and low energy minima 

If the minimum of V± does not lie close to a maximum of Vq there will be a smooth 
transition of the minimum to its final value. In B it was shown that for parabolic Vq 
and V\ this can lead to a substantial damping of the coherent modulus oscillations if C 2 
is large enough. 

In this section we study this damping in more detail for the gaugino condensation 
inspired potential with the minimum at T = 1.23 (and at the dual point T = 0.81). We 
assume the high energy potential V\ to be of parabolic form V\ = C 2 H 2 (T — T 1 ) 2 /2. 
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Fig 4-' T -trajectory for T\ = 1.001 and T\ = 1.7 for C = 10 



Unfortunately it is impossible to repeat the analytical approach of || to get the 
damping of the oscillation amplitude. This is, because the more realistic potentials we 
are interested in do not have the simple quadratic structure but contain the Dedekind 
function r] and its derivatives. 

Trying to approximate the model by an expansion which can be solved analytically 
proved to be inconsistent with the numerical solution because the terms one has to omit 
become too large. Therefore we had to stick to numerical simulations. 

We have studied the dependence of the final oscillation amplitude for a set of different 
starting values. The results are given in figs. 5 and 6. Each curve represents a number 
of runs with the same initial minimum T\ (given by the label at the curve) and shows 
the resulting (normalized) amplitude for different values of C on a logarithmic scale. The 
curve with the label 'Linde' represents the C-dependence of the oscillation amplitude of 
the model of ||. 

It is clear from these results that the damping of the moduli oscillations does not only 
depend on the value of C, but also strongly on the shape of the potential V between the 
minima of Vq and V\. If T\ is close to the final minimum, Vq can be approximated by 
a quadratic function. As expected, the reduction of the amplitude with respect to C is 
identical to the result in ||. However, outside of that region, there are strong differences. 
If Ti falls into a region of Vq which is rather flat, then increasing the value of C decreases 
the oscillation amplitude by a factor which is significantly smaller. If the modulus starts 
at a point where Vq is steep, then the reduction factor is larger. 
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Fig 5: Damping of oscillation amplitude in comparison with model of 
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Fig 6: Damping of oscillation amplitude for additional starting values 

These simulations consider C to be independent of T. However, as we showed in the 
last section, C will usually depend on T. It is therefore necessary to study the running of 
moduli with a non-constant C and one can expect that at least for T being in the vicinity 
of the low energy minimum, C should be large. 

The analysis confirms this: we choose T\ = 1.7 as the high energy minimum (starting 
value) and allowed C to vary quadratically from 1 to 8 between this initial and the final 
minimum like CiT) = 285(T - 1.7) 2 + 1 and C{T) = 285(T - 1.23) 2 + 1. Having C = 8 
at T = 1.7 and C = 1 at T = 1.23 (the second case) results in negligible damping of the 
modulus oscillation. The final amplitude is even larger than for constant C = 1 although 
over the whole trajectory of T the value of C is larger. 
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The first choice, where C becomes large as T — > 1.23, does show a damping. It 
corresponds however not to the maximum value of C along the trajectory but to some 
average. The final amplitude is comparable to the case of constant C ~ 2.5 for this 
example. 

To achieve significant damping it is therefore a necessary condition to have a large 
value of C around the final minimum T , and a consistently large C over the whole 
trajectory of the modulus T would be sufficient. 



3.4 Explicit examples for a smooth transition 

As we have shown in the preceding sections, a model with sufficient modulus damping 
necessarily will fulfill two criteria: first, the high energy minimum T\ must not be close to 
a maximum of Vo and second, C 2 (T ) = C 2 must be large. 

In addition it would simplify the analysis considerably if the trajectory of the modulus 
would be along the real axis. In principle this is not an unnatural requirement, since the 
theory obeys a symmetry T — > T and thus both Vq and V\ have an extremum on the real 
axis in the imaginary direction. This symmetry constitutes another example of how a 
symmetry of the theory could correlate the minima of the low and high energy potentials. 
Again the crucial point is to ensure that these extrema are indeed minima in the imaginary 
direction. 

Unfortunately within the specified class of theories we did not succeed in constructing 
such a 'real' model, which has large values of C 2 (T). However, models with T confined 
to the real axis and small C 2 (C 2 < 1) can be constructed. One example is F(T) = 
M^/( V "(T)j^(T)). 

Relaxing the condition for T to be real at all times we therefore attempt to construct 
models with large C 2 and complex trajectory of T. Furthermore we concentrate to on 
large final values Cq (the necessary condition), to get an impression of the size of C that 
can be expected. 

Using the Yukawa-couplings (|29"D one gets for C 2 at the two minima T = 0.81, 1.23 
(for broken supersymmetry) 

C 2 (1.23) = - J (-0.63 n + 4.49/ + 1.55) (44) 

q I. Ao — 1\ 

C 2 (0.81) = (-0.27/1- 6.78 1+ 1.33) . (45) 

q 0.81 — T\ 

One can see, that these models are in principle able to explain large values of C 2 , 
which is due to the properties of the modular functions (especially j(T)). The powers of 
the modular functions might however be constrained by the fact that V± should have a 
minimum at finite T. This is ensured if F(T) — > oo for Tr — > oo. This in turn gives the 
condition 

l>^-n-^. (46) 
12 12 y ! 
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One model which falls into this allowed region (modular weight n = —4) is then 

with C 2 (0.81) ~ 7/(0.81 - T x )/q. 

This is about as high as one can go if one stays within the orbifold limits on the 
modular weight (BD1) and 



4 Summary and conclusion 

Our paper discusses the behaviour of string moduli during the preheating period, the time 
after inflation when the energy density is dominated by the coherent oscillations of $. 

We show that $ can usually be integrated out to give an effective potential for the 
modulus, with a crucial dependence on the Yukawa coupling. 

For this paper we restrict ourselves to models which obey T-duality. As the low energy 
potential we use the one given by the standard gaugino condensation picture. For these 
models it is found that a stabilization of the modulus at high energy is possible for a 
limited range of the modular weight n of $ (n < 3) only. Within this limit, the reduction 
of the initial modulus amplitude around the high energy minimum can be calculated and 
will be 10 5 at most. To be compatible with standard nucleosynthesis one would need a 
reduction factor of about 10 10 . 

Studying the further time evolution of the modulus, we consider the different possi- 
bilities of how the modulus can settle to its low energy minimum. 

First we construct explicit examples with coinciding low and high energy minima at the 
self-dual points. In this case the initial stabilization will survive the subsequent evolution. 
Nevertheless the reduction is not large enough to completely solve the cosmological moduli 
problem, if the initial amplitude (at the end of inflation) is of the order 1 (in Planck 
units). One therefore needs a pre-stabilization by coinciding minima during inflation and 
preheating as well. 

Even then quantum fluctuations during inflation will drive the modulus away from its 
minimum, so that the additional reduction during preheating might be crucial to solve 
the moduli problem. 

A stabilization at self-dual points implies unbroken supersymmetry in the direction of 
the respective modulus. It can thus not work for all moduli, if supersymmetry is broken by 
a modulus at all. Another mechanism to avoid modulus oscillations, namely the smooth 
transition from high to low energy minimum (induced by a mass term C 2 H 2 with large 
C) is studied next. 

We found the modulus evolution in a more realistic (T-dual) potential to be more 
complicated than in a quadratic potential. If the high energy minimum is close to a low 
energy maximum the modulus oscillation is not significantly damped, irrespectively of the 
value of C 2 . This might be quite relevant, since the usual T-dual potential for broken 
supersymmetry has a saddle at the self-dual point T = 1. 
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The reduction factor will also depend on the shape of the low energy potential along 
the trajectory of T. If V is steep along this trajectory, the reduction will be larger for the 
same values of C. 

A nice feature of our framework is that it allows us to calculate C from the model. 
Large values of C are directly linked to steep couplings (large gradient with respect to 
T). 

Finally we analyze specific models and calculate the respective values of C(T). Our 
experience shows that it is difficult to construct examples which fulfill all requirements, 
namely: a well defined high energy minimum, a large value of C 2 (T) along the whole tra- 
jectory and (for convenience) a trajectory along the real axis. Though modular functions 
like j(T) are well suited to produce steep couplings, we could not obtain values of C 2 
larger than O(10). We believe that this in the upper region of what one can achieve. This 
can be contrasted to the values of C 2 needed to solve the moduli problem, which have to 
be at least one or two orders of magnitude larger. 

We believe that our work has serious implications for the cosmological moduli problem. 
The additional stabilization of the modulus at its high energy minimum during preheating 
will provide a large amplitude reduction which can help in solving the moduli problem. 
As we have shown, a further damping of the oscillation can be achieved via the two 
possibilities of coinciding minima and smooth transitions. For the second the oscillation 
amplitude will be governed by the factor C 2 {T) and the shape of the potential. However, 
it seems unlikely that the amplitude can be reduced down to the nucleosynthesis bound. 
Nevertheless, given some additional small reduction (e.g. by weak scale inflation), one 
might be able to achieve this. In any case these results can significantly lessen the demands 
on other mechanisms. 
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